Dynamics of matter-wave solitons in a ratchet potential 
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We study the dynamics of bright solitons formed in a Bose-Einstein condensate with attractive 
atomic interactions perturbed by a weak bichromatic optical lattice potential. The lattice depth is 
a biperiodic function of time with a zero mean, which realises a flashing ratchet for matter-wave 
solitons. The average velocity of a soliton and the directed soliton current induced by the ratchet 
depend on the number of atoms in the soliton. We employ this feature to study collisions between 
ratchet-driven solitons and find that soliton transport can be induced through their interactions. 
In the regime when matter-wave solitons are narrow compared to the lattice period the ratchet 
dynamics is well described by the effective Hamiltonian theory. 
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The ratchet effect, i.e. rectified average current in- 
duced by an asymmetric potential and unbiased zero- 
mean driving, has historically attracted a lot of atten- 
tion due to its possible relevance to biological transport, 
and molecular motors, and prospects for nanotechnol- 
ogy [11, i. Both in classical and quantum systems, the 
ratchet effect has been studied both in dissipative and 
Hamiltonian regimes [1, 01 and it appears to be due to 
a broken space-time symmetry of the perturbing poten- 
tial Directed transport can be experimentally im- 
plemented in a range of different physical systems rang- 
ing from semiconductor hetcro-structures to quantum 
dots, Josephson junctions, and cold atoms in optical lat- 
tices j6i]. Recently, the interest in optical ratchets and 
especially in the effect of interaction on the ratchet trans- 
port has resurged with the experimental advances 
in implement ing atomic ratchets for Bose-Einstein con- 
densates HI- 



As a physical system with intrinsically present non- 
linear interactions, Bose-Einstein condensates (BECs) of 
atoms with negative scattering length support the exis- 
tence of localized collective excitations - bright matter- 
wave solitons. It is therefore interesting to explore the 
effect of the ratchet potentials on the transport and in- 
teraction properties of such collective excitations. It is 
especially important because theoretical studies of the 
soliton ratchets so far were mostly focused on topological 
solitons whereas BECs represent a perfect test-bed 
for the study of the ratchet dynamics of a general class of 
non-topological solitons described by a continuous Gross- 
Pitaevskii (or nonlinear Schrodinger) equation. 

In this Letter, we study, for the first time to our 
knowledge, the effect of the ratchet potential on non- 
topological non-dissipative bright matter-wave solitons. 
A BEC represents an intrinsically nonlinear system which 
displays a single-particle behavior on a macroscopic scale 



being easily manipulated by reconfigurable optical and 
magnetic potentials. The ratchet potential is realized by 
means of a bichromatic optical lattice which is "flashed" 
on and off in such a way that its time-averaged amplitude 
vanishes. We show that both the ratchet effect and soli- 
ton directed transport are observed in this system even 
in the absence of losses, which sets it apart from previ- 
ously studied dissipative nonlinear systems subjected to 
ratchet potentials A weak potential does not af- 

fect the soliton shape in a wide range of parameters, and 
especially in the regime when the solitons are strongly 
localized, which justifies treatment of the soliton as an 
effective classical particle. However, we show that the 
ratchet still works even when the extended nature of the 
excitation cannot be ignored and the soliton width ex- 
ceeds the period of the lattice. Moreover, we investigate 
the influence of the ratchet on soliton scattering, and 
show that multiple collisions between solitons may pro- 
vide a space averaging mechanism that can enable obser- 
vation of a soliton current in a ratchet potential. 

To study the dynamics of bright solitons in attrac- 
tive BECs with ratchet potentials we consider a soliton 
formed in a strongly elongated condensate cloud [l4] sub- 
jected to a flashing one-dimensional optical lattice. As 
long as the energy of the longitudinal excitations is not 
sufficient to excite the transverse modes of the conden- 
sate, the system can be treated as one-dimensional and 
described by the Gross-Pitaevskii (GP) equation: 
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where the Fourier-synthesized lattice potential, 
V{x, t) ^ Vq f{t) [coa{x) + cos(2x + (f))] , 



(1) 



(2) 



is driven biperiodically: f{t) = sm{u!t) + sin{2ujt), and 
Vo depends on the intensity of the laser beams forming 



2 



the lattice. Quantum transport of atoms in the station- 
ary potential of the form ([2]) has been recently studied 
experimentally with ultracold rubidium vapor . With 
the time-dependent potential ^ the space inversion sym- 
metry is broken for (j) ^ 0,tt and the time inversion sym- 
metry is broken by our choice of f{t). According to the 
symmetry analysis [Hj , this can allow a directed transport 
of matter- wave solitons. 

In derivation of the dimensionlcss model ([T]) we have 
assumed that the energy, length, and frequency are mea- 
sured in the units of Eq = h'^k'^/m, uq = and 
ujo = hk'^ /m, respectively, where m is the mass of the 
atoms, and k is the wavevector of the optical lattice. We 
have furthermore reduced the three-dimensional mean- 
field model to the one-dimcnsional GP equation by as- 
suming that the wavefunction is separable in the follow- 
ing manner: ipsoix, y, z) = %l;iD{x)^{y, z), where $(y, z) 
is the ground state wavefunction, normalized to one, of 
the two-dimensional harmonic oscillator with the trans- 
verse trapping frequency LOi_. With these assumptions, 
the wavefunction ^E* in Eq. ^ relates to ijjio as fol- 
lows: vff = ipiD^giD, where gio = 2(asWj_)/(aoWo) is 
the renormalized interaction coefhcient that character- 
izes s-wave scattering of the condensate atoms with the 
scattering length a^. Number of atoms in the system is 
defined as: M — N/guj, where N — J \'i>\'^dx is the norm 
of the dimensionlcss onc-dimensional wavefunction. 

In the experimental setup [3] a bright soliton forms 
in the ^Li cloud with a modified scattering length 
Qs ~ —0.21 nm trapped in a quasi-one-dimensional 
atomic waveguide with uj± — 2tt x 710 Hz. We con- 
sider this experimental setup with an additional opti- 
cal lattice formed by light beams crossed at the an- 
gle = 38°, derived from a CO2 laser with the wave- 
length A = 10.62 //TO. With these assumptions the 
scaling units of length and frequency take the values 
Oo = A/[47rsin(6'/2)] = 2.52 /tm and tJo = Stt x 224 Hz. 
A stable bright soliton typically created in the experi- 
ment [l3| contains A/" ~ 5 x 10'^ atoms, which corresponds 
to iV R:! 2.62 in our units. We note that by changing the 
angle 9, it is possible to achieve smaller or larger values 
of ap, and hence of N, for the same value of A/". 

At t = 0, i.e. in the absence of the lattice, the shape 
of the soliton is given by the exact solution of the GP 
equation, 

xo, 0) = (7V/2)sech [(7V/2) {x - xq)] , (3) 

where xq is the initial position of the soliton, and the 
motion of the soliton is free. If we assume that the addi- 
tion of a weak ratchet potential does not affect the soliton 
shape during its evolution, then the Hamiltonian descrip- 
tion of the mean field allows us to treat the matter-wave 



the effective potential. 
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Using Eq. ^ and the equation of motion for the soliton 
position, cPxq/ dt^ = —dVcs/dx^, the soliton velocity can 
be obtained straightforwardly. The instantaneous shape 
of the effective potential, shown in Fig. [1] for f{t) = 1, 
becomes closer to the shape of the optical lattice, V{x, t) 
(Fig. [U thin solid line) , as N grows and the soliton be- 
comes more localized. For small A'' the second term in 
(|4]) becomes exponentially smaller than the first one, and 
hence the effective potential becomes practically symmet- 
ric, as demonstrated by a solid thick line in Fig. [TJ 




soliton as an effective classical particle [16| that moves in 



FIG. 1: (Color online) Effective potential Vcs vs. soliton 
position Xq, for = 1 (solid line), A = 3 (dashed) and 
A = 5 (dotted). Thin solid line is the optical lattice potential 
V{x,t). Parameters are: Vo — 0.3, (j) = n/2, f{t) — 1 

The effective-particle approximation (EPA) predicts 
that the cumulative velocity of the collective excitation, 
V = (1/r) Jp v(t)dt, is a function of the number of atoms 
contained in the soliton (A^). As can be seen in Fig. [21 
this prediction is in good agreement with the numerical 
results obtained from the GP model, where we have used 
T « 10'^ X 27r/cj. As Fig. [2] demonstrates, for a fixed 
initial position of the soliton, xq, there is a sharp transi- 
tion between a regime where a soliton oscillates between 
neighboring wells but is not transported [Fig. [2fb)] , and 
a regime where the soliton acquires a ballistic motion 
[Fig. [l^c)]. The simulations of ballistic motion are per- 
formed with the periodic boundary conditions, however 
we make sure that the tails of the soliton do not overlap. 
The EPA model also predicts that the velocity of the bal- 
listic motion tends to a constant value as A?^ — > 00 (and 
T 00), which is qualitatively confirmed by numerical 
simulations. 

The dependence of the cumulative velocity v on the 
number of atoms is also a function of the driving fre- 
quency, Lu [see inset in Fig. ([2])] . For high frequencies the 
velocity vanishes as w oc as the soliton becomes less 
and less affected by the rapidly oscillating potential. For 
the moderate driving frequencies used here the soliton 
does not radiate, and small variations of its shape do not 
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lead to particle loss. In the limit of small w the soliton 
strongly interacts with the lattice and may break up. In 
the effective-particle model the trajectory of the soliton 
in this case is chaotic and a detailed study of its behavior 
is beyond the scope of this paper. 




FIG. 2: (Color online) (a) Cumulative velocity v vs. num- 
ber of atoms in the soliton, TV, calculated using Eq.llJ (solid 
line) and EPA (dashed); v = 1 corresponds to 3.5 mm/s. Pa- 
rameters are: V = 0.3, — tt/2, uj — 10, xo ~ 0. Inset: 
Cumulative velocity v vs. driving frequency lj, for A'^ = 4 
and Xo = (solid line) and xq = —n/2 (dashed line). (b,c) 
Density plot of the mean field evolution, \il:{x,t)f, shown for 
the corresponding points at A = 2 and A = 4 in (a). 



The dependence of the soliton velocity on atom number 
is a general feature that occurs for symmetric periodic po- 
tentials as well, however the precise relationship between 
velocity and atom number depends on the initial position 
of the soliton relative to the lattice. Since the space in- 
version symmetry of the lattice is broken, we can expect 
that averaging over all initial soliton positions, xo, will 
lead to a directed soliton current. To demonstrate this 
effect, in Fig. |3]we plot the soliton velocity averaged over 
all initial positions, (v) = (l/27r) J^^ vdxo, as a function 
of the number of atoms. Although the average velocity is 
always non-zero, one can identify two different regimes of 
the ratchet dynamics depending on the number of atoms, 
as discussed below. 
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FIG. 3: (a) Average velocity {v) vs. number of atoms in the 
soliton. A, calculated using the GP model (solid line) and 
EPA (dashed). Parameters are: Vo = 0.3, <j) = tv /2, cu = 10. 



As seen in Fig. |3l for small values of A'' the EPA re- 
sults and numerical soluton of Eq.([T]) disagree both on 
the onset of the ratchet effect and on its magnitude. For 
A^ < 2.5 the soliton's size is comparable to or larger than 
a period of the optical lattice. Hence it is more accu- 
rately described as a wavepacket than an effective par- 
ticle. The details of the soliton response to the flashing 
potential are best seen by examining the dependence of 
the cumulative velocity v on the initial position of the 
center-of-mass, xq. In Figs. [4l^a,b) we show this response 
for A^ = 1 and N = 2. Interestingly the numerical results 
show that the soliton has either no cumulative velocity, 
or moves in only one direction. As a result, for A^ < 2.5 a 
soliton attains a much larger average velocity than that 
predicted by the EPA theory (see Fig. [3]). In fact the EPA 
incorrectly predicts that the soliton can move in both di- 
rections depending on its initial position. We also note 
that for A^ = 1 the cumulative velocity expected from the 
EPA is almost symmetric around xo = as expected in 
the limit of small A^ due to the symmetry of the effective 
potential (see Fig[T]). 
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FIG. 4: (Color online) (a) Cumulative velocity, v, vs. soliton 
initial position, xo, for (a) A = 1; (b) A = 2; (c) A = 5; 
calculated from the numerical solution of the GP equation 
(solid line) and EPA (dashed), (d) Poincare section for the 
eff'ective-particle model in (c). Parameters are: Vo — 0.3, 
(j> = 7r/2, LO = 10. 

For large values of A^, corresponding to a strongly lo- 
calized matter-wave soliton, the ratchet dynamics is well 
described by the EPA. In Fig. [3] we observe a good agree- 
ment between the ratchet velocity predicted by the effec- 
tive particle model and that found from solving Eq.([T]) 
numerically. Similarly, in Fig. IH^c) we see a good agree- 
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ment between the numerics and the EPA in the de- 
tails of the cumulative velocity dependence on xq. In 
Fig. md) we show the Poincare section corresponding to 
the effective particle results of Fig. UJc) . Two different 
types of trajectories are observed, transporting or non- 
transporting. A comparison between Figs. and d) 
shows a clear correlation between cumulative velocity in 
Fig. m^c) and either non-transporting trajectories (corre- 
lating with the zero- velocity) or transporting trajectories 
(correlating with either positive or negative velocity). 

As shown above, solitons with different N travel at 
different ratchet-induced speeds even if their relative ini- 
tial position in the lattice is the same. Multiple soliton 
collisions could therefore be realized if the ratchet poten- 
tial was combined with a toroidal trap [l3|- As can be 
seen in Fig. El^a) a larger moving soliton can then induce 
the transport of a smaller soliton which would otherwise 
not be transported by the ratchet. This is due to the fact 
that each collision incrementally changes the soliton's po- 
sition in the phase space and its original motion can thus 
be dramatically modified if it is eventually moved from a 
non-transporting to a transporting trajectory. The driv- 
ing has little effect on the actual near-elastic soliton col- 
lisions because of its small amplitude and fast variation. 

If the solitons have equal values of N, collisions can 
occur only if they have different relative positions in the 
lattice, and therefore different speeds. In this case how- 
ever, the interaction between the solitons is strong due 
to the large interaction energy. Each collision therefore 
induces a much more pronounced transition to a differ- 
ent phase-space trajectory. This is evident in Fig. El^b) 
where we can see dramatic changes in the soliton veloc- 
ities after collisions. In this scenario, the spatial shift 
that solitons acquire during each collision may lead to an 
effective averaging over xq after multiple collisions. As a 
result, a nonzero total average current can, in principle, 
be observed for a sufficiently large number of collisions 
or for a sufficiently large number of interacting solitons. 

In conclusion, we have studied the ratchet dynamics of 
bright solitons in Bose-Einstein condensates and demon- 
strated that solitons may either be transported through 
the potential or oscillate, depending on their initial po- 
sition relative to the lattice. For small atom numbers, 
the soliton transport occurs in one direction only, while 
larger solitons may be transported in either direction. 
We show that the rate of transport for a given initial po- 
sition is atom- number dependent, with solitons that con- 
tains more atoms typically moving faster. Importantly, 
we find that averaging over all soliton initial positions re- 
veals an overall directed soliton current, and show a good 
qualitative agreement between the numerical results and 
an effective particle approximation. Finally, we illus- 
trate how the ratchet potential affects interactions be- 
tween solitons containing the same or different number 
of atoms. Collisions can cause an instantaneous tran- 
sition between different trajectories in the phase space. 
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FIG. 5: (Color online) (a) Collision between two soliton, one 
with A*' = 4 with initial position xq = and the other one with 
= 2.2 and initial position xq = 47r. (b) Collision between 
two soliton with = 4 located at a;o = and xo = Stt -I- 1.2. 
Parameters are: Vb = 0.3. (j> — it/2, lu — 10. 



including that between non-transporting and transport- 
ing trajectories. This effect could potentially be used for 
directed transport or spatial filtering of solitons based on 
the atom number. This work is supported by the Aus- 
tralian Research Council (ARC). 
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